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. Using effective field tlieory, we investigate nuclear modification of nucleon parton distributions 
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Abstract 



(for example, the EMC effect). We show that the universality of the shape distortion in nuclear 
parton distributions (the factorisation of the Bjorken x and atomic number {A) dependence) is 
' model independent and emerges naturally in effective field theory. We then extend our analysis to 



. study the analogous nuclear modifications in isospin and spin dependent parton distributions and 

o , 

\ generalised parton distributions. 
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In 1983, the European Muon Collaboration (EMC) reported measurements Q| of the 
ratio Rpeix) of F2{x) structure functions in iron and deuterium in deep inelastic scat- 
tering (DIS). In the parton model, these structure functions are defined as -F^(x) = 
YliQl^ [ifi^) where qf{x)lqf{x)] is the parton distribution function (PDF) for 

quarks [ant i- quarks] of flavor z in a nucleus. A, and Qi is the electric charge of g^. Bjorken x 
is the longitudinal momentum fraction of the struck parton in the infinite momentum frame 
with respect to its parent nucleon. Rpf.{x) is normalised to unity if the nucleons in the 
nuclei are non- interacting. The main observation of EMC was that Rpei^) deviated from 
unity by up to 20% over the range 0.05 < x < 0.65 in which measurements were taken. The 
EMC result was unlooked for and came as a surprise to many physicists at the time partly 
because the typical binding energy per nucleon is so much smaller (<1%) than the nucleon 
mass and the energy transfer involved in a DIS process. 

Over the past two decades, further experiments have been performed by many groups 
(see Refs. j^lsl 0, 0] for recent reviews) aiming to better understand the details of nuclear 
modifications of hadron structure functions (here referred to generically as the EMC effect). 
The EMC result has been confirmed and demonstrated in many other nuclei ranging from 
helium to lead. Nuclear modification of structure functions has also been studied in other 
situations such as proton-nucleus Drell-Yan experiments and quarkonium production [s,!^. 
Overall, a very interesting picture has emerged; for an isoscalar nucleus of atomic number 
y4, the shape of the deviation from unity of Ra{x) = F^{x) / AF^ (x) (here has been 
converted to the isoscalar subject to a small model dependent error, and the slight 
dependence of Ra{x) is suppressed) is universal 0, ^, namely independent of A within 
experimental error bars, while the magnitude of the distortion is empirically proportional 
to the number density of the nucleus, pa |9| . Fits to the available data that support these 
features will be presented elsewhere [lO|. These findings have inspired many theoretical 
analyses seeking to understand the details of the EMC effect in various approaches (see the 
reviews for summaries). Different physical processes have been identified as the causes of 
the modifications in different x regions; e.g., nuclear shadowing at low x and Fermi motion 
at large x. A recent conclusion of detailed model studies was that the EMC effect necessar- 
ily implies modification of the nucleon PDFs and cannot be explained through traditional 
nuclear physics [11]. However, little is known directly from QCD. 

In this Letter, we employ effective field theory (EFT) to investigate the EMC effect by 



2 



studying nuclear matrix elements of the twist-two operators which are related to parton 
distributions and structure functions via the operator product expansion. We find that the 
universality of the shape distortion of the EMC effect is a model independent result, arising 
from the symmetries of QCD and the separation of the relevant scales. The x dependence of 
Ra{x) is governed by short distance physics, while the overall magnitude (the A dependence) 
of the EMC effect is governed by long distance matrix elements calculable using traditional 
nuclear physics. We then proceed to study analogous nuclear effects in isospin-odd, and 
spin-dependent parton distributions and generahsed parton distributions (GPDs) [l^. We 
also discuss aspects of extracting the shape of the EMC effect from first principles using 
lattice QCD. This approach provides a clear connection between the EMC effect and many 
other observations of nuclear modification of hadron properties. 

EFT is a model independent approach which only makes use of the symmetries and 
scale separation of the system (for recent reviews see Ref. ^^). This approach has been 
successfully applied to many low energy processes in A = 1,2,3,4 systems. Recently EFT 
has been applied to the computation of hadronic matrix elements of twist-two operators 
in the meson and single nucleon sectors 14 , Q| and app lied to chiral extrapolations of 



lattice calculations of moments of parton distributions jl_7|. The approach has also been 
extended to analyse moments of generalised parton distributions jl8!], large Nc relations 
among PDFs in nucleons and the A-isobar 15|, and deeply virtual Compton scattering 
in the nucleon |3| and deuteron systems 1. The method is readily generalised to the 
multi-nucleon case. Although we will concentrate on quark bilinear twist-two operators, 
the framework can be easily applied to gluonic operators and thereby to nuclear effects in 
gluonic distributions which are important in heavy ion collisions at RHIC and LHC and at 
a future Electron-Ion Collider. 

To described the EMC effect observed in F2 data on isoscalar nuclei, we consider the 
normalised, spin singlet, isoscalar twist-two operators, 

0M-Mn ^ g^(Mo^/)A.i . . . iD^'-)q/ (2M"+i) , (1) 



where (...) indicates that enclosed indices have been symmetrised and made traceless, D 
{D^ — D^)/2 is the covariant derivative and M is the nucleon mass. The matrix elements 
of in an unpolarised single nucleon state with momentum P can be parametrised as 

^p|(^MO...Mn|p^ _ ■ ■ ■i^'^", (2) 



where the nucleon velocity = P^/M. It is well known that the coefficients {x"')q corre- 
spond to moments of the isoscalar combination of parton distribution functions, 

{x^)^ = l^xx-qix), (3) 

where q{x) is the isoscalar quark distribution and q{—x) = —q{x). 

We first consider only nucleonic degrees of freedom {i.e., assume that pions are integrated 
out of the EFT - they will be reintroduced below) and perform the standard matching proce- 
dure in EFT, equating the quark level twist-two operators to the most general combinations 
of hadronic operators of the same symmetries 
hadronic operators in the matching are 



;wo operators ; 



The leading one- and two-body 



Q^O-f^n ^ (a;")^^;MO . . . ^;Mn^t^ + UnN^ N] + ■ ■ ■ , (4) 

where v'^ = v'^ + 0{1/M) is the velocity of the nucleus. Operators involving additional 
derivatives are suppressed by powers of M in the EFT power-counting. In Eq. we have 
only kept the SU(4) (spin and isospin) singlet two-body operator anv'^° ■ ■ ■ v^" {N^N) m the 
above equation. The other independent two-body operator f3nV^° ■ ■ ■ v^" {NWN) , which IS 
non-singlet in SU(4) [r is an isospin matrix], is neglected because /5„/a;„ = 0{1/N^) ~ 0.1 
21 1, where is the number of colors. Furthermore, the matrix element of (^NWnY for 



an isoscalar state with atomic number A is smaller than that of (A^^A^)^ by a factor A 
. Three- and higher- body operators also appear in Eq. numerical evidence from 
other EFT calculations indicates that these contributions are generally much smaller than 



two-body ones 
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Nuclear matrix elements of O^o-- /^" give the moments of the isoscalar nuclear parton 
distributions, qA{x). The leading order (LO) and the next-to-leading order (NLO) contri- 
butions to these matrix elements are shown in Fig. 1(a) and 1(b), respectively. For an 
unpolarized, isoscalar nucleus, 

= v^"" ■■■v^'"{A\0^'-^'"\A) (5) 
= {x^),[A+{A\an{N^Nf\A)], 

where we have used (y4|A^'''A^|y4) = A. Notice that if there were no EMC effect, the w ould 
vanish for all n. Also aQ = because of charge conservation. Asymptotic relations [23[ and 



analysis of experimental data 



24| suggests that cti ~ 0, implying that quarks carry very 
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FIG. 1: Contributions to nuclear matrix elements. The dark square represents the various operators 
in Eq. Q and the light shaded ellipse corresponds to the nucleus, A. The dots in the lower part 
of the diagram indicate the spectator nucleons. 

similar fractions of a nucleoli's and a nucleus's momentum though no symmetry guarantees 
this. From Eq. © we see that the ratio 

rv 

(6) 
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Or, 



is independent of A which has powerful consequences. In all generality, the isoscalar nuclear 
quark distribution can be written as 

qA{x) = A[q{x)+g{x,A)]. (7) 

Taking moments of Eq. Eq. (jUj) then demands that the x dependence and A dependence 
of g factorise, 

g{x,A)=g{x)g{A), (8) 

with 

g{A) = {A\{N^Ny\A)/AAl, (9) 



and g{x) satisfying 



1 



dx x"'g{x) 



(10) 



Aq is an arbitrary dimensionful parameter and will be chosen as Aq = 1 fm~^. Crossing 
symmetry dictates that the even and odd a„ separately determine the nuclear modifications 
of valence and total quark distributions. These results apply to any isoscalar combination 
of parton distributions including ^2(0;) for isoscalar nuclei. Thus our result implies that 



Ra{x) = = 1 + g^^ix)giA), (11) 

which says that the EMC effect (the deviation of Ra{x) from unity) has an universal shape 
described by gF2{x) while the magnitude of the deviation, G^A), only depends on A. 

The above analysis gives a simple explanation of the observed universal shape of the EMC 
effect, or equivalently, the factorisation ofg{x, A). The key to establishing this factorisation 
is that other sources of nuclear modification contributing to the right-hand side of Eq. Q 
must be suppressed (higher order in the EFT) such that the A independence of Eq. © 
can be established. We stress that the factorisation persists when pions are included in our 
analysis. In Fig. ^ examples of the leading pionic contributions are shown. The various 
single-nucleon diagrams, such as Fig. ^c), simply renormalise the nucleon moments, {x"')q, 
without contributing to the EMC effect. Two- and more- nucleon diagrams such as those in 
Fig. Hfd) andme) contribute to the EMC effect, but only at N^LO and higher (see Ref. Q 
: or explicit calculations) since the pions must be radiation pions rather than potential pions 
2^ (we are free to choose the twist-two indices to be /ij = for all i). Other contributions 
that could upset the factorization include a two-body operator which is similar to that in 
Eq. but with two more derivatives. However this operator also contributes at N^LO. 
Consequently, the universality of the EMC effect is preserved to good accuracy. For large 
X it is clear that the factorisation must break down (simply consider the region a; > 2 in 
which only three- and higher- body operators contribute) though the structure function is 
very small in this region anyway. We stress that the factorisation is a model independent 
result and is just a consequence of scale separation and the SU(4) spin-isospin symmetry 
from large A*";,. 

It is clear from Eq. Q that Q{A) is governed by long distance physics which can be 
computed using a traditional, non-relativistic nuclear physics approach. It is interesting to 
note that the mean field scaling ofQ{A) ^ pA ^ log(^) describes the empirical A dependence 
of Ra{x) well , even though the mean field approximation is not justified for nuclei 

where the two-particle S-wave scattering lengths are much larger than the mean distance 
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between nucleons. 

Information on the shape distortion function g{x) is encoded in the short distance pa- 
rameters a„ associated with the strength of the two-body currents. One can either fix the 
an from experimental data (to determine all a„, data on F,^{x) and F^{x) are required) 
or calculate {NN\0^'^'"^'^\NN) in two nucleon systems to extract them. The latter ap- 
proach, however, is intrinsically non-perturbative and thus requires lattice QCD. Recent 
analysis j^l suggest that one can use background fields coupling to twist-two operators to 
extract information on the a„ from the finite volume scaling of two particle energy levels 
0,123, At present, only the first few coefficients could be extracted on the lattice 

because of problems with operator mixing 3Cj. However, even a calculation of these would 
be significant since it would be addressing nuclear modification of hadron structure from 
first principles. Lattice calculations can also be used to investigate the large Nc prediction 
of /5„/a„ = 0{1/N^) discussed below Eq. OH). 

Given the success of the EFT approach in explaining aspects of the isoscalar, helicity 
averaged EMC effect, we shall now proceed to study the isospin and spin dependent cases. 
For the isovector operators 

(^Mo-M„ _ q^i^^o^D^^l . . . iD^'-'^T^q/ (2M"+^) , (12) 

EFT operator matching leads to 

(^Mo...A.„ = {x''),X° ■ ■ ■ v^-N^T;N [I + as.niV^AT] + ■ ■ ■ , (13) 

where the ellipsis includes the higher order nucleonic and pionic operators. Analysis similar 
to the isoscalar case implies 

" {Z -"^^)t(x) = ^ + ^'^^^ ' 
where qz{x) = u{x) — d{x) and Z{N) is the proton (neutron) number. To test this fac- 
torisation, one can either consider the difference between F2's in {Z,N) = {n + m,n) and 
(n, n + m) mirror nuclei and compare it with F2 — F2, or disentangle ua{x) and dA^x) 
with the proposed neutrino- nucleus experiment, MINERi/A 3^ . 

To generalise our analysis to spin dependent operators is also straightforward. The op- 
erators related to moments of quark helicity and transversity distributions are 

^AgJ" = qi^^^l^iD^"' ■ ■■iB^^\^ql (2M"+^) , (15) 
O^^o^-Mn _ g(jP{M^5^^Mi . . . iD^'-h^q/ (2M"+2) . 



In the isoscalar case (tq = 1), the matching yields 

^AgJ" = 2(x")AgA^"^5('^''w^^ ■ ■ ■ w^-^AT [1 + 7„,oA^^A^] + ■ ■ ■ , (16) 

and similarly for the other operators (the spin operator S'^ = (0, cr/2) where cr are Pauli 
spin matrices). Again, we have similar factorisation; e.g., 

where AA^ is the differences between the number of nucleons with positive and negative spin 
projections in the longitudinal direction. Similarly, Qa{A) = (A|A^^A^A^^cr3A^|A)/AA^AQ. 
Whilst there is data on longitudinal asymmetries in light nuclei from which the gi{x) struc- 
ture functions can be extracted, disentangling nuclear effects in the unpolarised and po- 
larised structure functions will be difficult. Currently, polarised heavy nuclei targets, in 



which modifications would be larger, are not available. Recent model calculations |22| find 
nuclear effects in the polarised structure function gi{x) to be significant. Nothing is known 
experimentally about the transversity structure function even in the proton but analogous 
nuclear modifications can be derived. 

It is also possible to study nuclear effects in GPDs by computing off-forward matrix 
elements of twist-two operators jl^. In the quark contribution to nuclear spin, JgA, for 
example, in addition to extending the operators matched to 0^°'^^ in Eq. by replacing 
V — > v+i-^ (reparameterisation invariance Q] constrains this form), we also need to consider 
the term 

-2:^^D^^^N[S'^'^,S,](^v + ^^Y N[l+vN^N]^ , (18) 
where Jq(g)N is the quark(gluon) angular momentum content of the nucleon, to obtain 

J,A = {x)qNLz + 2.J,^Sz [1 + 7] ni^A)] . (19) 

Remarkably, explicit calculation shows that the Lz term is free from two-body current 
corrections even though ai 7^ in Eq. (j3)). Similarly, the gluon contribution to nuclear spin 
satisfies 

JgA = {x)gNLz + 2Jg^Sz [l-v n{A)] , (20) 

where the same constant rj appears by total angular momentum conservation. Consequently, 
using the sum rules {x)qN + {x)gN = 1 and JgN + JgN = 1/2, we recover 

JqA + JgA = Lz + Sz . (21) 
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For details, see 




To summarise, we have studied the EMC effect (nuclear modification of parton distribu- 
tions) in EFT and seen that the scale separation of short and long distance effects provides 
a model independent derivation of the factorisation of the x and A dependencies of Ra{x), 
relying only on the symmetries of QCD. Similar factorisations are predicted to occur in other 
probes of nuclear structure such as spin-dependent structure functions and GPDs. 
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